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Abstract. LEP and SLC provide accurate data on the process e e~ — > // at the Z resonance. 
The GigaZ option at a future linear e + e~ collider (ILC) will further improve these measurements. 
As a consequence, theory predictions with sufficiently smaller errors are necessary in order to fully 
exploit the experimental accuracies and to derive indirect bounds on the scales of new physics. 
Here we review the currently most accurate predictions of the Z pole observables (e.g. sin 2 O^g, rz, 
a had) m the context of the Minimal Supersymmetric Standard Model (MSSM). These predictions 
contain the complete one-loop results including the full complex phase dependence, all available 
MSSM two-loop corrections, as well as all relevant Standard Model contributions. 

PACS. f2.60.Jv Supersymmetric models - 12.15.Lk Electroweak radiative corrections 



1 Introduction 

Z boson physics is well established as a cornerstone of 
the Standard Model (SM) [1]. Many (pseudo-) observ- 
ables [2] have been measured with high accuracy at 
LEP and SLC using the processes (mediated at lowest 
order by photon and Z boson exchange) 



a particularly challenging task) . As a step in this direc- 
tion we compute the currently most precise predictions 
for the observables at the Z resonance [3] . These con- 
tain the full one-loop result, all available higher order 
MSSM terms and all relevant SM contributions. For 
the first time we include the full phase dependence at 
the one-loop level. 
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at a center of mass energy y/s ~ Mz- In particular 
these arc the effective leptonic weak mixing angle at 
the Z boson resonance, sin 2 9 e s, Z boson decay widths 
to SM fermions, F(Z — > //), the invisible width, il nv , 
the total width, Fz, forward-backward and left-right 
asymmetries, Afb and Alr, and the total hadronic 
cross section, cr° ad . Here we focus on sin 2 9 e g and rz, 
as these two observables show the strongest sensitivity 
on effects due to virtual SUSY particles [3]. 

Together with the measurement of the mass of the 
W boson, Mw, and the mass of the top quark, m t , the 
Z pole observables have been instrumental in bound- 
ing the mass of the SM Higgs boson, the last free pa- 
rameter of the model. In a combined fit containing 
Z pole observables, W mass and W decay width, the 
indirect constraints predict a SM Higgs boson mass 
of Mjj = 76±g GeV, with an upper limit of Mu < 
144 GeV at the 95% C.L. [1]. 

In order to fully exploit the high-precision measure- 
ments, the theoretical uncertainty in the predictions of 
the (pseudo-) observables should be sufficiently smaller 
than the experimental errors (in view of the antici- 
pated ILC precisions [4, 5] for the Z observables this is 
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2 The Z pole observables 

Z pole pseudo observables are commonly defined and 
calculated in an effective coupling approach. This ap- 
proach exploits the fact that the dominant contribu- 
tions to the process e + e~ — > // at s ~ M§ stem from 
resonant Z boson exchange diagrams. Non-resonant 
terms arise from photon exchange diagrams and box 
contributions, both of which are accounted for as part 
of the unfolding procedure from the experimental data 
(see Refs. [2,3] for details). The electroweak radiative 
corrections can thus be absorbed into effective vector 
couplings, g v , and axial vector couplings, g A . These 
effective couplings are in turn used to define effective 
fermionic mixing angles (at Born level these coincide 
with the weak mixing angle, sin(? w = s w ) 
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The latter are commonly expressed as [2] 

ry = N£r \ Pf \ (4(7/ - \K f \) 2 R f v + R{) , 



with 



1 — 
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and derived from the axial vector couplings as de- 
tailed in Ref. [3]. In the above equations, Qf , I3 and 
N[ stand for the charge, third isospin component and 
colour factor of the respective fermion /. The radiation 
factors Ry A account for QED and QCD interaction of 

the final state fermions in the decay Z — ► //. Tau lep- 
ton and more importantly bottom quark mass effects 
also enter via R VA - The effective mixing angles defined 
in eq. (2) are intimately realted to the forward back- 
ward and left right asymmetries measured at LEP and 
SLC. The total Z boson width is obtained by summing 
over all partial widths ry 
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3 One-loop result and higher order terms 

The computation of our MSSM predictions for the 
observables at the Z resonance consists of two main 
steps: the computation of the full MSSM one-loop re- 
sults, for the first time under consideration of CP- 
violating complex MSSM parameters, and the inclu- 
sion of all available higher order contributions from 
SM and MSSM. Details of the calculations, which are 
briefly summarised in the following, can again be found 
in Ref. [3]. 

To accomplish our first goal we renormalise the 
Zff vertex in the relevant parameters and compute 
all contributing MSSM one-loop graphs and countert- 
crms. All relevant Feynman graphs are calculated mak- 
ing use of the packages FeynArts [7] and FormCalc [8]. 
As regularisation scheme dimensional reduction [9] is 
used, which allows a mathematically consistent treat- 
ment of UV-divergences in supersymmetric theories 
at the one-loop level. The inclusion of loop corrected 
Higgs masses and couplings in the complex MSSM is 
another new feature of our calculation. Our implemen- 
tation is performed in accordance with the program 
FeynHiggs [10,11] and the discussion in Ref. [12]. We 
furthermore resum the tan enhanced bottom Yukawa 
couplings following an effective coupling approach [13]. 
For the first time we perform a full one-loop calcula- 
tion for decay of the Z boson into ncutralino pairs, 
Z — ► XiXi) which contributes to the invisible width 
and consequently also the total width of the Z boson, 
provided m^o < M^/2. 

In Ref. [3] we give an exact description of the higher 
order terms which are included in our predictions for 
the observables at the Z resonance. Our philosophy re- 
garding the inclusion of higher order terms is strongly 
influenced by fact that the theoretical evaluation of 



the Z pole observables in the SM is significantly more 
advanced than in the MSSM (see Ref. [14] for a re- 
cent discussion of the state-of-the-art results in the 
SM). In order to obtain the most accurate predictions 
within the MSSM it is therefore desireable to take 
all known SM corrections into account. This can be 
done by writing the MSSM prediction for a quantity 

x = 9v,A>Pl> K f>- ■ ■ as 
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where x SM is the prediction in the SM with the SM 
Higgs boson mass set to the lightest MSSM Higgs bo- 
son mass, M Hl , and x mssm-sm = x susy d cnotes the 

difference between the MSSM and the SM prediction. 
In order to obtain x MSSM according to eq. (7) we eval- 
uate 2; MSSM - SM a t the level of precision of the known 
MSSM corrections, while for x SM we use the currently 
most advanced result in the SM including all known 
higher-order corrections. As a consequence, x takes 
into account higher-order contributions which are only 
known for SM particles in the loop, but not for their 
superpartners (e.g. two-loop electroweak corrections to 
Ak beyond the leading Yukawa contributions). Besides 
including all known higher order SM contributions, we 
also account for all available generic SUSY two-loop 
terms [15,16], which enter as universal corrections via 
the p parameter. 



4 Dependence on complex parameters 

As already mentioned, sin 2 c r and Fz are the two 
Z observables which show the strongest sensitivity on 
effects of new physics. We therefore focus on these 
two observables in the following. Analyses including 
a® ad , Ri, Rb etc. can be found in Refs. [3,6]. Detailed 
numerical studies showed [3] that the dependence on 
the sfermion mass parameters is much stronger than on 
the chargino/higgsino parameters. We therefore only 
investigate the dependence on the phases of A t and 
A b . As for ArjM w [17], we find that the effective cou- 
plings glf^ only depend on the absolute values \X t \, 

\Xf,\ of the off-diagonal entries in the t and b mass ma- 
trices, where X t — A t — \ij tan/3, X b = — \i tan/3. 
Thus, the phases of p, A t and Ab only enter in the 
combinations (4>A t b + </v)) giving rise to modifications 
of the squark masses and mixing angles. It furthermore 
follows that the impact of 4>A t {<pA b ) on the sfermion 
masses is stronger for low (high) tan 0. 

In Figs. 1 we show sin 6 e & and Fz as a function 
of (j>A t (with (j>^ — (pA b — 0), for different values of 
tan/3 (varied from tan/3 = 5 to tan/3 = 45). Shown 
in the green-shaded bands are the current experimen- 
tal values in their la range. Using the same conven- 
tions as in Refs. [3,17], the other parameters are set to 



Ms = M H ± = M 2 =m- g = 500 GeV, \A tA , 



1000 GeV, 4> Ml = 4>m 2 = 4>g = 0. 

As expected, the dependence of sin 2 9 e Q and Fz 
on (f>A t is most pronounced for small tan 0. The vari- 
ation of <pA t in this case gives rise to a shift in the 
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two precision observables by l-2a. The effect becomes 
smaller for increasing tan/3, up to tan/3 =15. On the 
other hand, for high tan (3 the lighter b mass becomes 
rather small for the parameters chosen in Figs. 1, reach- 
ing values as low as about 100 GeV for tan/3 = 45. 
This leads to a sizable shift of ~ 1— 2a in the Z ob- 
servables already for vanishing phases. The slight rise 
in the dependence on cf>A t for tan /3 > 25 is due to 
the overall enlarged SUSY contributions which occur 
for large tan [3 and the resulting low sbottom masses. 
We checked that the dependence on (f>A t is in general 
significantly larger than the dependence on 4>A h ■ 
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Fig. 1. Prediction for sin 9 e g and Fz as function of the 
phase of the trilinear coupling At- The other SUSY param- 
eters are: Mi = M H ± = M 2 = = 500GeV, A T = A t — 
A b = fj, = 1000 GeV, </) M = <j)M 1 = <j>M 2 =4>g = 0, </>A b = 



5 MSSM parameter scans 

We now investigate the behaviour of sin 2 9 e g, the Z ob- 
servable most sensitive to higher order corrections, by 
scanning over a broad range of the SUSY parameter 
space. The SUSY parameters given in Tab. 1 are varied 
independently of each other, within the given range, 
in a random parameter scan. Unlike Refs. [6,18,19], 
where our results for the electroweak precision observ- 
ables were already employed in CMSSM multiparam- 
eter analyses, the scans in this section are entirely un- 
constrained, apart from the fact that we require the 
Higgs mass bounds and direct SUSY particle exclu- 
sions limits from LEP searches to hold. 

The SM and the MSSM predictions for sin 2 9 cS as 
a function of m t , obtained from the scatter data with 
m t as an additional free parameter, are compared in 





Parameter 


Range 


sleptons 


Ml, , M E , 


100 to 2000 GeV 


squarks 


M Qz , M Vi , M Dz 


100 to 2000 GeV 




At, A b 


-2000 to 2000 GeV 


gauginos 


Mi,M 2 


100 to 2000 GeV 




M 3 


195 to 1500 GeV 


Higgs 


M A 


90 to 1000 GeV 




tan/3 


1.1 to 60 






-2000 to 2000 GeV 



Table 1. MSSM parameter ranges in the unconstrained 
random parameter scan shown in Figs. 2 and 3. Family 
indices i run from 1 to 3, where the respective parameters 
are varied independently. All phases are set to zero. 

Fig. 2. The predictions within the two models give rise 
to two bands in the mt-sin 2 9 e g plane with only a 
relatively small overlap region (indicated by a dark- 
shaded (blue) area). The allowed parameter region in 
the SM (the medium-shaded (red) and dark-shaded 
(blue) bands) arises from varying the only free parame- 
ter of the model, the mass of the SM Higgs boson, from 
M| M = 114 GeV, the LEP exclusion bound [20] (lower 
edge of the dark-shaded (blue) area) , to 400 GeV (up- 
per edge of the medium-shaded (red) area). The very 
light-shaded (green), the light shaded (green) and the 
dark-shaded (blue) areas indicate allowed regions for 
the unconstrained MSSM. In the very light-shaded re- 
gion at least one of the ratios m^/m^ or fn-^Jm^ 
exceeds 2.5 (we work in the convention that m; < 
m 7 ) , while the decoupling limit with SUSY masses 
of 0(2) TeV yields the upper edge of the dark-shaded 
(blue) area. Thus, the overlap region between the pre- 
dictions of the two models corresponds in the SM to 
the region where the Higgs boson is light, i.e., in the 
MSSM allowed region (M h < 130 GeV [10,11]). In 
the MSSM it corresponds to the case where all super- 
partners are heavy, i.e., the decoupling region of the 
MSSM. The 68% C.L. experimental results for m t = 
(170.9± 1.8) GeV [21] and sin 2 6 cS are indicated in the 
plot. As can be seen from Fig. 2, the current experi- 
mental 68% C.L. region for m t and sin 2 9 c s is in good 
agreement with both models and does not indicate a 
preference for either of the two. The prospective ac- 
curacies for the Tevatron/LHC (5 sin 2 ^^ vatron / LHC _ 
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2 fl ILC/GigaZ 



5m 



ILC/GigaZ 



0.1 GeV) are also shown in the plot (us- 



ing the current central values), indicating the strong 
potential for a significant improvement of the sensitiv- 
ity of the electroweak precision tests [5] (see Rcf. [22] 
for a recent review of the anticipated errors at future 
colliders). 

In Fig. 3 we show the combination of My/ [17] and 
sin 2 9 c s with the top quark mass varied in the range of 
165 GeV to 175 GeV. The ranges of the other varied 
parameters and the colour coding are the same as in 
Fig. 2. The current 68% C.L. experimental results for 
My/ and sin 2 9 e g are indicated in the plot. The region 
of the SM prediction inside todays 68% C.L. ellipse 
corresponds to relatively large m t values, outside the 
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current experimental range of mt — (170. 9± 1.8) GeV. 
Thus, the combination of Myy and sin 2 9 c s exhibits 
a slight preference for the MSSM over the SM. The 
anticipated future improvements in the measurements 
of M\y and sin 2 9 e s are again indicated. 
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Fig. 2. Unconstrained MSSM random parameter scan for 
sin 2 # e ff as a function of mt over the ranges given in Tab. 1. 
Todays 68% C.L. level ellipses as well as future precisions, 
drawn around todays central value, are indicated in the 
plot. 
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Fig. 3. Unconstrained MSSM random parameter scan over 
the ranges given in Tab. 1 and mt = 165 . . . 175 GeV. 
Shown is the combination of Mw and sin 2 f9 c ff. Todays 
68% C.L. level ellipses as well as future precisions, drawn 
around todays central value, are indicated in the plot. 
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